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ABSTRilGT 


Utilization of the earth’s magnetic field for 
attitude control of gravity gradient satellites is studied. 
A controller consisting of a single hody-fixed onboard 
electromagnet is proposed for achieving simultaneous 
control of the three-dimensional rotational motion of the 
spacecraft. The Liapunov stability theorem is applied to 
the linearized system equations to synthesize a general 
control law for th'e dipole strength. A response analysis 
establishes the validity of the approach in orbits with 
different inclinations. The system offers enhanced 
reliability as no moving pants are involved. The seai- 
passive character of the S7 promises an increased 
operational life-span for the satellite. 
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1. INTaOIUGTIOH 

The general motion of a space vehicle is quite complex. 
It may he thought of as the superposition of two distinct 
modes of motion: the translation of the vehicle centre of 
mass and the rotation of the spacecraft about the centre of 
mass. The curvilinear translation of the mass centre relative 
to the primairy attracting body describes the orbit of the 
satellite. The rotational motion about the mass centre, on 
the other hand, determines the orientation of the vehicle in 
space. It is also referred to as the attitude motion of 
the satellite. 

The success of a vast majoriiy of space applications, 
such as communications, weather, military, earth-resources, 
etc., depends on the ability of the spacecraft to maintain a 
preferred orientation in space. Unfortunately, a satellite 
initially oriented along the desired attitude deviates from 
it due to several reasons. These are: the gravity gradient 
effect, magnetic field interactions, aerodynamic forces, 
solar radiation pressure, uncompensated motion of internal 
machinery and micrometeorite impacts. Magnetic field 
interactions constitute a significant perturbing effect for 
satellites below 1600 tan altitude. Solar radiation pressure 
represents the principal disturbing effect for hi^ altitude 
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missions, e.g. , in the geostationary orhit. The effect of 
aerodynamic forces is relatively unimportant for satellite 
altitudes above approximately 300 km. 

Por successful orientation of the satellite, an 
attitude control system is therefore required. The various 
methods available for attitude control may be classified as 
active, passive and semipassive according to the power 
consumption involved. The use of control elements sudi as 
gas jets, reaction wheels and gyros come under active methods 
of attitude control. The active methods require a large 
energy consumption v/hich limits the life-span of the satellite. 
A high speed of response and accuray can be obtained by these 
techniques. 

Passive methods, on the other hand, do not involve 
any energy consumption. Gravity gradient stabilization and 
spin stabilization fall in this category. In gravily gradient 
stabilization, the values for the moments of inertia about the 
principal axes of the satellite are so chosen that the vehicle 
will have gravitational restoring torques about all the three 
axes. The vehicle then tends to orient itself with its axis 
of minimum moment of inertia along the local vertical and 
the maximum moment of inertia axis normal to the orbit plane. 
Spin stabilization technique is based on the fact that a 
spinning body exhibits gyroscopic stiffness against external 
torques. The satellite is spun about its axis of symmetry 
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and it tends to maintain its angular momentum vector fixed 
in inertial space* The passive methods, in general, lead to 
an increased operational lifetime of the satellite as no 
energy consumption is required. Their pointing accuracy 
capabilities, however, are limited. 

The semipassive methods represent a compromise 
between the active and passive methods. The concept involves 
utilization of suitable controllers whictt interact with the 
space environment to generate control torques about the sate- 
llite centre of mass. The utilization of the aerodynamic 
forces, solar radiation pressure and magnetic field inter- 
actions, come under this class. In aerodynamic control, 
suitably designed control surfaces are extended from the 
satellite to intercept the relative atmospheric flow to 
obtain the control torques. Mirror-like surfaces may be 
used at high altitudes to generate control moments from solar 
radiation pressure. In magneti c attitude control system, 
electromagnetE are used onboard to interact with the earth’s 
magnetic field. "By varying the onboard dipole strength and 
orientation, one can govern the magnetic torques about the 
centre of mass of the satellite. The energy consumption in 
these methods is limited to activating the controller 
elements. The semipassive methods provide a high pointing 
accuracy and a moderate speed of response. They neither 
involve any mass- expul si on nor the use of eleaents requiring 
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a large power consumption, and hence, promise an increased 
life-span for the satellite. 

Generation of control torques through the interaction 
of onboard electromagnetic dipoles and the earth’s magnetic 
field appears to be particularly attractive as the system 
reliability is enhanced by the elimination of moving parts. 
This thesis deals with some aspects of the magnetic attitude 
control of spacecraft. Important previous investigations in 
the area are briefly discussed next. 

Brief Survey of Previous Work: 

Libra tion damping of gravily -gradient spacecraft 

using passive hysteresis damping was x)nsidered by Alper 
1 2 

and O’Neill and Pischell . The authors suggested using 
rods of high permeability magnetic alloys for damping the 
attitude motions. The response rates of the hysteresis 
dampers, however, are extrmely slow. The concept of dosed- 
loop magnetic torquing and its feasibility was suggested by 

■x 4 

White et al. Banium and Mackison applied the sample and 

hold concept to a gravity gradient satellite employing three 

mutually perpendicular onboard electromagnets for effecting 

5 

three-axis attitude damping. Wodi and Pande investigated 
the feasibility of three axis nutation damping and attitude 
control of both spinning and nonspinning s.-- tellites using 
a hybrid magnetic-solar control system. The authors 
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coasidered two magnetic controller models involving a 
rotatable electromagnet and two fixed electromagnets* for 
spin axis orientation, and the use of solar radiation 
pressure for simultaneous control of the satellite pitch 
attitude. 

The problem of maintaining the spin-axis of a 
spinning satellite normal to the orbit plane has been a 
subject of considerable discussion. Ergin and Wheeler^ 
proposed an autonomous magnetic attitude control system for 
the purpose employing a dipole along the pitch axis. 

7 

Subsequently, Wheeler presented a feedback control law 
in oluding provisions for active magnetic nutation damping 
as well as directing the spin axis of an axisymmetric 
spinning satellite in a circular orbit toward any direction 
in sjace. The stability of the control law was demonstrated 
via the method of Liapunov. However, a torque-free environment 

s 

was assumed. Renard obtained an open loop control law that 
requires quarter orbit switching for near polar orbits. 

Q 

Rlugge Lots, and Maltz considered a contactor control system 

» 

to stabilize the attitude of a spinning space vehicle and 

obtained the system response using a linear switching 

criterion. The authors accounted for gravity torques during 

the analysis but the system stability was not demonstrated. 

Sonnabend”*^ and I^hl and Poxman^^ developed precession control 

12 

laws that utilize a roll axis dipole. Tossman proposed a 
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contx’ol ^stem for magnetio nutation and despin control wliicli 
utilizes the concept of enhanced magneti c hysteresis damping. 

-T -j X 

Tratblik et al . ’ proposed an attitude control system for the 

spin-axis orientation of the lES-4 and IjES-5 spacecraft using 

a system of electromagnets as the sole source of corrective 

15 

torque at the synchronous altitude. Sorenson developed 

magnetic pointing systsn ap|)li cable to spacecraft with a 

1 6 

limited attitude detemnination capability. Shigehara 

developed a control law based on the asymptotic stability 

criterion for the spin axis and spin rate control using 

dipoles along the spin axis and perpendicular to it, 

17 

respectively. Collins and Bonello discussed several 

magnetic control schemes that provide nutation as well as 

1 8 

precession control. Alfriend proposed a closed loop 

control law for a dual spin satellite using a dipole 

19 

parallel to the spin axis. Alfriend and Stickler developed 
a three axis closed loop attitude control system for Earth 
observatory momentum bias spacecraft for initial acquisition, 
precession control, nutation damping and pitch axis momentum 
control. 


It is apparent that the possibility of achieving 
simultaneous control of the roll, yaw and pitch motions by 
means of a single, body -fixed electromagnet has not been 
investigated so far. On the other hand, it would represent 
a significant reduction of the hardware design without any 



7 


loss of system reliability. This thesis proposes and 
develops a magnetic attitude control system for gravily 
gradient satellites involving a single electromagnet with 
a body -fixed orientation in the spacecraft principal axis 
system. A control law for the dipole strength, promising 
asymptotically stable operation, is ^nthesized and the 
performance of the control system is examined throu^ a 
response analysis of the system. 
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2. POMULATIOR OF THE PROBLM 

2.1 Geometiy of Motion and Controller Configuration: 

Fig, 2. la shows the orbital motion of an un symmetrical 
satellite with its centre of mass S moving in a Keplerian 
orbit around the centre of the earth 0. Pg and P^ indicate 
the geographi c and magnetic axes of the earth directed from 
the earth's centre to its geographic and magnetic TTorth 
poles, respectively. The two are offset by an angle e 
( - 11,4°). As the earth spins about its geographic axis, 
the intersection of the line of intersection of the plane 
containing P„ and P„ with the equatorial plane advances 
uniformly with time. 

The rotational motion of the satellite about its 
centre of mass is shown in Fig. 2. 1b. To describe the 
satellite attitude relative to the rotating orbital 
co-ordinate system x^, the following sequence of 

finite rotations may be employed: 

(i) the rotation B about the x^ axis bringing x^^, y^, 
to x>j , y-^ » z^ , 

(ii) the rotation yabout the y^ axis bringing x^jy-jtZ^ 

to X2>y2*^2* finally, 

(iii) the rotation X about the Z 2 axis bringing 

in coincidence with the body -fixed satellite principal 


axes x,y,z. 
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The rotations , and are referred to as the jaw, 
roll and pitch attitude angles, respectively. Any spatial 
orientation of the spacecraft nay he uniquely described by 
allowing appropriate yaw, roll and pitch rotations in their 
proper sequence. 

Pig. 2,2 shows a schematic representation of the 
proposed magnetic controller which consists of a single 
electromagnet with its dipole moment P, Ihe mag»itude of F depend 

upon the area enclosed by the current carrying coil and the curre: 
passing through it. It# direction is determined by the 
direction of the current in the coil according to the 
right-handed screw rule. The unit vector p indicating the 
normal to the plane of coil is located by its spherical 
Co-ordinates 5^ and ^ relative to the principal axes x,y,z . 

2.2 Equations of Motion: 

The equations of motion for the spacecraft may be 
derived using either the Lagrangian or the Eulerian 
formulation. The overall system possesses five degrees 
of freedom, namely, R, 6 , 3 , Y and X , The generalized 
co-ordinates R, 9 describe the orbital motion while 0 ,y 
and X describe the attitude motion of the satellite. The 
two, in general, are coupled and this leads to a very 
complicated five degree-of-freedom systen. However, it is 
well Known that the influence of the attitude motion on 
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the orbital motion is negligible^ Therefore, one may 
directly derive the equations of attitude motion and 
simplify the derivation significantly making use of the 
fact that the orbit remains essentially Keplerian. 
Presently the derivation for the case of a circular 
satellite orbit v/ith angular speeds = flis given using 
the Eulerian approach. 

2.2.1 Eulerian formulation: 


Euler’ s equations of motion for the rotational motion 
of a rigid body about its centre of mass are veil known to 
b e 


0} 

+ 

dz-Iy) 

0) 

CO 

= 

(2.1a) 

X 


y 

z 

X 

(Jj 

+ 

(I -T-J 

0) 

CO 

= K 

(2.1b) 

y 


^ X z^ 

z 

X 

7 




0) 

CO 

= 

(2.1c) 

z 


X 

y 

z 


where the moments of inertia, angular velocity components 
and the external torque components are referred to the 
principal axes 3:.y,z. 


The external torque components consist _ 
of the contributions from the gravity gradient torque and 
the magnetic torque, ' i.e. , 


M = M + M 
“x gx mx 


M = M + M 

V gy ffiy 


{2.2a) 


(2.2b) 
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+ M 


mz 


(2.2c) 


The gravitational torque about the centre of mass arises 
due to the variation of the gravity field across the dim 
sions of the satellite. Its components are given hy^^ 


M 


= 3(-^) (ly-l 2 )sinY cos ysin X 


gx 


(2.3a) 


M = 3(-^-^)(l^-Ij2)siny cosycosX (2.3h) 

R 

M = 3(-^)(I -I )sin>. cosX cos^y (2.3c) 

where terms of th e order (satellite dimension/R)^ and 
higher are ignored compared to unity in its expansion. 

The magnetic torque generated by the interaction 
of the onboard electromagnet with the earth’s magnetic 
field may be written as 

M^=PxB=upxB (2.4) 

where 

u = I A (2.5) 


with I representing the current and A the total area 
enclosed by the current carrying coil. The unit vector 
p along the normal to the plane of the coil is simply 

P = Pj, I + Py 3 + Pz E 


C2.6) 



■where 


% 


= cos 

6 ^ cos 62 

(2.7a) 

= sin 

5-] cos ^2 

( 2 . 7 h) 

= sin 

2 

( 2 . 7 c) 


The geomagnetic induction vector B may he expressed in 
terms of its components in the x^, y^, system hy 

21 

assuming the earth to he a canted magnetic dipole as 


= (•|)sini sin2Ticose + (■|)[cos0 + (^)C (.1+cosi)cos( 2n”0) 
+ (1 -cosi) cos (2 t)+0)> ] sine (2.8a) 

B /Dj^ = (■|)sini(1-3cos2p)cose + (~)[cosi sin0 
yxi 

+ (|){ (lH.cosi)sin{2Ti-0)+(1-cosi)sin(2Tl+52!)}]sin£ 

(2.8h) 

B /I) = -cosi cose + sini sine sin0 {2.8c) 

zn' m 

where the geomagnetic dipole moment is 

D = M /R^ (2.9) 

in 6' 


with 

2^5 3 

M = 8.1 X 10 Gauss-cam . 
e 

Fete that the angle 0 varies due to the earlh's rotation 
according to 0 = n^. With 6 as the independent variable 
in the circular orhit, this gives 
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0 = (llg/ q) + 0^ 


( 2 . 10 ) 


The present analysis ignores any nodal regression as a 
dynamic effect which is equivalent to assuming the earth 
to he gravitationally spherical. If desired, it may he 
accounted for hy appropriately updating 0^ after every few 
orhits. The body components of B may he obtained from the 
components along the y^, systeaa through a series of 
co-ordinate transformations. The final expression is 
obtained as 




r 

cosX cosycosT] 

cosX cosysint) 

sinXsinp ^ 

Bx 

= 

-sinX cosPsinT] 

+sin XcospcosTi 

-sinycosX . 



-sinycosX sinpsinT) +sinYCOs Xsinp cost) cosp 

; 



-sinX cosycosri 

-sinXcosysinq 

cosXsinP 

' 

= 

-cosXcospsinp 

+C0S ^cospcosT) 

+sinXsiny. ' 



+ sin^sin^siny . 
sinn 

-sinXsinysinp . 

COST] 

cosp ; 


= 

sinycosTi 

+cosysinpsinT] 

sinysinn 

» cosysinpcosTi 

cospcosy 





(2.11) 


The magnetic torque is thus obtained in terms of its 
components in the principal axis system by substituting 
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EqrLS.(2.1l) and (2.7) in Eqn.(2.4) and carrying out the 
vector cross product. Subsequent substitution of these 
components into Eqns. (2.2) finally leads to the desired 
total external torque expressions. 

The determination of the kinematical relationships 
for ■‘-be body components of the angular velocity is fairly 
straightfor-ward. Again, the use of a number of co-ordinate 
transformations is involved, finally leading to 

0 ) = (y -i-^ sinp )sinX + (pcosy - Qcosp sinY)co^ (2.12a) 

0 ) = (y sinp)cosX - ( p cosy ^ cospsiny) sin X (2.12b) 

• • 

oi ^ =X + p siny +!^ cosp cosy (2.12c) 

Substituting Eqns. (2.1-2) in Eqns. (2.1) and introducing 
the total torque components, the equations of motion governing 
the attitude motion result. Eor convenience, the orbital- 
angle 9 = f^t may be introduced as the independent variable. 

The nondimensionalized equations of motion are finally 
obtainned as 
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p cosy - [k;^{ ^sinycosY sinx+ (X’ y’ + X*sir^) cos X 

+ X'sinX (cospsiny - p’cosy) + (p'y* + p'sinp)siny cos X 
+ P'sinx ( cosp sin^y — p’siny cosy) 

+ (y * + sinp)cosp cosy cos X + (cos^p cosy siny 
“ p‘cosp cos^y) sinX} cos X + K^{ 3sinY cosy cosX 
“ (x'y* + X 'sinp )sin X + X' cosX ( cospsiny-p ’ cosy ) 

- (P’y* -i- P 'sinp )sinysinX + p * cosX( cospsin^y 

- P'siny cosy) - (y ’+siiip )cosp cosy sinX 

2 2 
+ (cos p cosysiny - p'cosp cos y)cosX} sin X 

- P’sinp siny + P'y’ siny + y’cosP cosy -X *sinp-X*y ’] 



(2.13a) 



3sinY cosy sinX+(X ’ y'+^' siiip)cosX + 

+ x'siii>(a «BsinY - P’ cosy) + (P’y' + P’ s±np)sinY cosx 

+ P’ sinx(cosp sin^Y-P' sinY cosy) + (Y’+siJip )cosp cosycosx 

P 2 

+ (cos PcosY sinY-p ’ cosp cos Y)siiiX} sinX 

+ 3sinY COSY cos X- (X ’Y’ + ^’sinp)sinX 
4-X’cosx(cosp sinY-p’ cosy) (p 'y* +P ’ sinp )sinY sinX 
+ P’CosX(cosp sin^Y “ P' sinY cosy )~(y ’ )cosp cosy sinX 
+ (cos^p COSY sinY - P’cosP cos^y)oosX} cosX 
+ P* cosp + X'cosp sinY - P’y’oosY 


A 

{(p B - p 'B„)u sinX} 


+ — { (p^B^ - P^2^ ^ (2.13B) 


I 


x" 4 - ) [ 3sin^osXcos^Y + (Y’+sinp)^ sinXcosX 

'' 1.+E-J ^2 

- (Y*+sinp)(cospsinY - p ’ cosy) cos 2A 

- (cosp sinY - P' cosy)^ sinX cosX ] - p' sinp cosy 

- Y'cosp sinY + P'y' cosy + P" siny 


I 




{2.13c) 
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where the parameters J ^ and K2 defined as 





(2.1.4a) 



(2.1 4h) 
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3. CONTROL SYNTHESIS 


The objective of this chapter is to arrive at a 
suitable control strategy for the magnet! c moment of the 
onboard electromagnet. The general problem of establishing 
a Control policy for the control variable u is rather 
formidable in view of the complexity of the governing 
system of equations. An approach based on small amplitude 
motion near the equilibrium configuration of the system is 
therefore adopted. 

3.1 System Equilibrium and Stability: 


In absence of the control u, it is observed that 
3 = Y= 0 represents an identical solution of Eqns. (2.13). 
Physically, it implies that in absence of any external dis- 
turbances the principal axes x,y,z remain aligned with the 
orbital co-ordinates x^, y^, z^, respectively. This 
corresponds to an earth-pointing spacecraft attitude. IXi e 
to the effect of disturbances such as might result from 
initial misalignment or due to micrometeorite impacts, the 
attitude would deviate from its desired configuration, 
linearizing the equations of motion about the equilibrium 
point, they simplify to: 


3 " - 6 - (1+K^) Y 



(p B - p B )u 
'^t'y z y 


1 


(3.1a) 
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y" + 422 Y- (Kg-DB^ = (5.1b) 

7 

X" + 3(E;^+E;2)/(1+E^K2^= jJT~^Px®y “ (5.1c) 

Z 

TTobe that the components of the earth’s magnetic induction 
vector B^ are retained as sudi because these can be 

directly sensed by means of three mutually perpendicular 
magnetometers aboard the vehicle. 


The conditions for the stability of the equilibrium 
of an uncontrolled satellite can be obtained easily frcm 
Eqns. (3.1) by letting u = 0 and examining the characteristic 
equation 

/ 22 3(K^+Kp) 

{ s^ + (1+3Kp-K.K^)s -4K.£p}{s‘' + i —} = 0 (3.2) 

(UK^K^) 

The equation appears in the factorized form as the pitch 
motion of an uncontrolled satellite decouples from the 
roll-ya-w motion near the equilibrium position. It is 
apparent that one can at best expect stability in the sense 
of Liapunov corresponding to purely imaginaiy roots of the 
characteristic equation. Also, because of ,the constraints 
on the irfcertia properties of a rigid body given by 


I + I > I , I, 

X y z’ y z ^ z X ^ 

the parameters and £2 limited to the 



(5.3) 
PUR 


J. ^ ' 


m 


58344 
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-1<E^<1, -1<K2<1, £^£^>0 (3.4) 

In view of these inequalities, the range of , Kg over 

which stable equilibrium exists can be expressed in the 
22 

K-j-Kg space as shown in fig. 3.1. Normally, the moments 
of inertia of a gravily gradient satellite are designed 
such that the systan configuration in the K^, Kg space is 
located within the stable range shown. In the following, 
the satellite mass distribution is assumed to correspond to 
a nominally stable configuration. 

3.2 Control Law Derivation: 

The objective of the magnetic controller is to provide 
attitude control torques which maintain the earth pointing 
attitude of the satellite. A suitable control law, developing 
appropriate torques in all the three degrees of freedom, is 
presently synthesized via the second method of Liapunov 
>ased on the linearized equations of motion. The construction 
of a suitable Liapunov function becomes relatively easy by 
expressing th e coupled roll-yaw equations in the state-variable 
form. 


Letting x^ = B , Xg = P' ? ^3 " Y » ^4 = ‘Y’ » 

Ejuations (3.1a and b) can be expressed in the state- 


variable form as 






24 



( 5 . 5 ) 

A transformation matrix which ■uncouples the roll-y'^'w motion? 
of ah uncontrolled satellite can he found using methods of linear 
systan- theory. One such transformation is described by 


y = Q X (5 .6) 



( 3 . 7 ) 


Here, the characteristic frequencies roll -yaw 

motion may be obtained as the solution of ttie equation 
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4 

0 ) 


(1 + 3K2 - K^K2) J' - 4K^£2 = ° 


(3.B) 


When the transformation described by Eqns.(3.6) is applied 
to Eqn.(3.5), only the coupling in the terms involving the 
con tool variable v/ould remain. Explicitly the result is 
found to be 




(3.10a) 


(3.10b) 
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A Liapunov function for the complete system is now easily 
selected as 

V = yf + yj + “(y, + y|) + + (5^^) (J-n ) 

where a > 0 is an arbitrary constant. 


Differentiating Eqn.(3.11) gives 

1 + 

v = 2[ y^y^l + ^2^2 + x‘ + (3(1^ +£—)>'" 1 

( 3 . 12 ) 

Substituting for the derivatives from Eqns* (3*9) (3-10), 

all terms not involving the control variable cancel out 
leading to 

O 

of"] 4-OCG^ ) 

y, = L- u[(A-,p + 


I. 


“x 


(PyB^ - P^By) + AjCp^y - PyB,,) X ■] (3- 13 ) 

where the constants A-j, A 2 » A^, A^ are defined by 

q 1 . ,°= (“f ^ 

'^l " ^^l+ac2 '^”1 +Ki 


A, 


I 

■y 

1.0 


, 4K,aC^ 

i w__A — A + -A- ) 


“2 ■ (1-K,) 


(3.1-4a) 

(3-14b) 

(3.t4c) 


^ = (V)(T^2n“ ♦ 

^ -^y 1 +aC2 


(5*14d) 
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1 

1+a c| 


1 + 

)( ^ 


) 


(3.l4e) 


Eqn.(3*'i3) 2 Eiy te utilized to determine a control law for 
u which, promi.qes asymptotic stability of the system by 
ensuring tha.t V’ < 0. 


Eow an error function a is defined as 
<J = (A - PA) + (V + V)(Vz ■ 

+ - PyB^) (5- ■'5) 

and asymptotic stability is assured if 

u=f(a) (3-16) 

where f(cf) is any function vrith the property 

of(0) < 0 (3*17) 

Eqns. (3*1-6) and (3.1. 7) determine the control law for the 
onboard dipole moment. Both linear as well as nonlinear 
control laws follow from these. The linear form is written 
as 

u = -E o, E > 0 (3*18) 

where E is a constant. 

A typical nonlinear form would be the bang-bang control 
law 

u = -Uq sgn 0 (3-19) 














where repre^^'ents the maximum amplitude of the magnetic 
moment at which the controller operates all the time. 

The mechanization of the proposed control system 
would involve the use of three magneto-meters oriented 
along three mutually perpendicular axes and attitude error 
as well as rate sensors. Besides, appropriate function 
generators would he required to continuously determine 
the control signals. A control flow diagram depicting the 
control system is shown in Fig. 3* 2. 
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4. RESULOB AND DISCUSSION 


In order to evaluate the performance of the control 
syster, the full nonlinear equations of motion (2.13) were 
solved in conjunction with the control law synthesized using 
the linearized system. As pointed out earlier, an analytical 
solution of the equations is not possible and therefore numerical 
integration was undertaken. A digital computer program using 
the Hamming’s predictor-corrector method with a Runge-Kutta 
starter was developed. When the linear control law was 
Considered, a step size of 0.5® gave results of sufficient 
acoiracy over a wide range of qystem parameters. 


As the system behaviour is characterized by a large 
number of xarameters, it would take an execessive amount of 
computational effort to investigate the effect of each of 
them. Therefore, a specific satellite geometry was selected 
for the numerical study. The parameter values taken are: 


I^ = 80 kg-m 


I = 120 kg-m 


= 150 kg-m 


2 




-0.3750, Eg = 0-5333 


Orbit inclination i = 45° 

Orbit altitude R = 1000 km. 
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The associated control law was chosen to le the linear 
form (eqn. 3.18) and the controller parameters taken as: 

<5^ = 45°, §2 = ^ 

These lead to the constants (i = 1,5, 4, 5) occurring ia the 
error function to he 


= 0.0079, A^ = -0.0740 

A^ = 0.9007, A^ = 0.5873 

The response of the ^stem subject to initial position 
and velocity errors was investigated first. These correspond 
to attitude errors which might result from an initial 
pointing misalignment and micrometeoritic impacts, respec- 
tively. The ^stem resx>onse to a simultaneous initial error 
of 5° in roll, yaw and pitch is indicated in Pig. 4*1. ^he 
time-histories of the control variable u and the decay of 
the Liapunov function V are also shown. It is observed 
that all the three degrees of freedom are practically 
damped out (< 0.1°) in about 5 orbits. The maximum onboard 
magnetic moment commanded during the process is found to be 
1 .58 amp-m^. Similar results when the initial errors are 
doubled to 10° are presented in Pig. 4. 2. Here, all the 
three degrees of freedom settle within 0.1° in about 7 orbits. 
A comparison of the two cases suggests that the maximum 
magnetic moment required is almost double of the previous case 






ii# 
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On the other hand, the system damping rate is higher ■when 
the initial errors are larger. This is also reflected in 
the nature of the decay rate of the V function -whidi is a 
measure of the overall systaa damping. In this sense the 
behaviour of the present nonautonomous system appears to 
he similar to that of a linear time invariant system. 

figures 4.3 and 4.4 show the system response for 
the case of a correctly positioned satellite subject to 
impulsive disturbances. Again, simultaneous disturbances 
are considered in all the three degrees of freedom. The 
responses clearly indicate the capabilily of a single 
magnetic torques in countering the dis'fcurbance. It not 
only removes the resulting errors but the maximum attitude 
excursions during the process are also not too large. The 
associated dipole strength required is also quite snail. The 
system settling times to within 0.1° are found to be approxi- 
mately 5 and 12 orbits for the two initial conditions 
considered. It may be pointed out here that the disturbance 
values of 0.1 and 0.2 represent an extremely severe situation 
deliberately chosen to evaluate the controller performance 
under the most adverse conditions. In reality, the impulses 
due to micrometeorite impacts are at least an order of 
magnitude lower than these. Consequently, the performance 
of the system would be expected to be far superior to the 
responses presented here. 




K» 






Fl9»4.4 SysteiR response and control history in an inclined orbit 
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It -was decided next to examine the effectiveness of 
the approach in orbits with different inclinations. The 
results for the case of i = 45 ° have already been discussed. 
Typical responses for the case of a polar orbit (i = ^>0®) are 
shown in Pigs. 4. .5 and 4.6 for initial position and velocity 
disturbances, respectively. The damping rates appear to be 
somewhat lower than for the case of i = 45°. The effect of i 
on the behaviour, however, is not easy to isolate since the 
electromagnet locations 6 -j , 62 and the inclination i together 
determine the nominal angle between the. dipole axis and the 
geomagnetic induction vector B at a given orbital position. 

As the latter determines the distribution of the total control 
torque between the three degrees of freedom, the performance 
in any orbit would improve through an optimum choice of 
and <5 2. An exception occurs for orbit very dose to the 
equatorial plane (i= 0 ) where the pitch component of the 
control torque remains small for any 6 .|, ^2 because the 
vector B is nominally at e = 11.4° from the pitch axis. 
such a case, the overall effectiveness of the control system 
reduces as indicated by the responses shown in the Pigs. 4-7 and 
and 4.B. The system performance was found to be quite 
sensitive to the choice of K and a. It appears possible to 
reduce the settling times significantly through an optimum 
selection of these parameters. The nature of the equations 
being nonautonomous, an optimization study would have to 
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consider the response due to distirhances that might occur 
at any position in the orbit. The effort required appears 
Vorthvfhile only for a specific configuration during the 
design stage. 

Finally, a comment about the effect of the orbit 

I 

altitude would be important. When the Eqns. (2.3, 2.9 and 

2.10) are substituted in Eqns. (2.12), the orbit radius R 

2 3 

is absorbed in the constant 1 /q R = 1/p . The magnetic 
torque components then contain R explicitly, only in the 
terms of order e which is small. Therefore, the non- 
dimen si onalized responses over a range of R- would remain 
practically the same except when i “ 0. 

An attempt was also made to evaluat e the response 

of the system in conjunction with the bang-bang control law 

2 

(Eqn. 3.19). With u^ = 5 amp-m , a = 1.5 and other 
parameters as before, considerably faster damping was 
indicated by the numerical results. Unfortunately, sufficient 
quantitative agreement could not be found between the various 
computer results obtained with step sizes ranging from 
1° to 0.01°. This appears to be due to the high sensitivity 
of the qystem to the switching times of the control u. Efforts 
to improve the quantitative agreement were found to involve 
a prohibitive amount of computation time. 
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5. CO^'I.JLUDI^G SSCmXS AlH) RSCOMIIETOi'^TIOITS 
FOR FUTURE WORK 

The main conclusions based on the study may be 
summarized as follows: 

1 . A magnetic controller employing a single electro- 
magnet has been proposed which is capable of damping the 
roll, yaw and pitch disturbances simultaneously. 

2. ii feedback control law for the onboard electro- 
masJnet has been derived using Liapunov’ s second method of 
stability. The control obtained is of a general nature 
pemitting a wide choice of the dependence of the control 
variable on the error function. 

3. The ^stem is found to be qute effective against 
both initial attitude errors as well as impulsive disturbances. 
Its performance may be further improved throu^ an optimum 
choice of the controller parameters. The systen appears suitable 
for applications not demanding a high speed of response, such as, 
long-life scientific satellites. 

4. The ability of the controller in orbits with 

any inclination has been demonstrated, althou^ the response 
is relatively slower in near-equatorial orbits. 

5. The system does not involve any mass expulsion 
schemes. Only a small amount of electrical energy is required 
to activate the electromagnet which can be generated onboard 
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ty means of solar cells. As there are no moving parts, it 
represents an inherently more reliable approach, il’h'as the 
controller promises an increased life-span for the spacecraft. 

The following recommendations are made for future work 
based on the present study: 

1. It would be desirable to obtain more precisely 
the results with the proposed bang-bang control law. This 
would also permit an optimization of the system performance 
with respect to the weighting parameter 

2. The procedure should be extended for the case of 
a satellite in an elliptic orbit. At high altitudes during 
an elliptic orbit, the earth's magnetic field weakens 
further while the solar radiation pressure assumes a 
dominant role. It appears interesting to develop a hybrid 
magnetic- solar controller permitting attitude control by a 
semi-passive approach throughout the orbit. 

3. At low altitudes, both the aerodynamic forces 
and the earth's magnetic field are significant. A magnetic 
control approach could be developed to counter the continuous 
disturbance resulting from the aerodynamic torques. 
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